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C^ . Abstract 

We report on a scaling test of several mesonic observables in the non-perturbatively 
0(a) improved Wilson theory with two flavors of dynamical quarks. The observables are 
constructed in a fixed volume of 2.4 fm x (1.8 fm)^ with Schrodinger functional boundary 
conditions. No significant scaling violations are found. Using the kaon mass determined 
in [1], we update our estimate of the Lambda parameter to A^/mx = 0.52(6). 



1 Introduction 

In this article we summarize the results of a set of simulations of QCD with two de- 
generate flavors of quarks employing Schrodinger functional boundary conditions [2]. 
The range of quark masses covered corresponds to a ratio of the pseudoscalar mass to 
the vector mass, Mpg /My, in the interval [0.4, 0.75]. Our final goal is to compute the 
fundamental parameters of perturbative QCD, namely the scale parameter A and the 
quark masses Mq, in units of a hadronic observable such as the Kaon decay constant 
-Fk- We emphasize our effort to control all systematics. Here we focus on cutoff effects 
and reach (for one quark mass) a lattice spacing that is smaller than those previously 
achieved in large- volume simulations of the 0(a) improved Wilson action [HElllllS]. 

While simulations of QCD with at least A^f = 2+1 flavors of sea quarks are mandatory 
to provide accurate non-perturbative predictions with direct phenomenological implica- 
tions, in our view the N{ = 2 theory represents a framework well suited to address a 
number of fundamental aspects of low-energy QCD that have not been clarified yet, a 
couple of which we shall presently mention. 

One such question is the A^f dependence of A-^ / Fk and Mg/ Fk- Since these quan- 
tities have been computed in the quenched theory [SI [7j , it is interesting to know the 
separate effects of the (up, down) quarks and those of the strange quark. To our knowl- 
edge, the influence of the strange sea quarks on hadronic observables has not been 
demonstrated very clearly so far. 

Secondly, it is important to determine the quark mass at which one-loop SU(2) chiral 
perturbation theory becomes accurate at the (say) 3% level. We see a strong motivation 
to address this question in the N{ = 2 theory, with one parameter less to tune on the 
QCD side. And with a small number of low-energy constants in the chiral perturbation 
theory, this is probably the cleanest way to establish the latter as the low-energy de- 
scription of QCD from first principles. Given the level of accuracy one is interested in, 
all sources of systematic error have to be addressed. In particular any observed non- 
linearity in the quark-mass dependence of Fpg and Mpg must first be shown to survive 
the infinite volume limit before it can be claimed that the chiral logarithms have been 
observed. Cutoff effects represent an additional source of systematic uncertainty, which 
is computationally expensive to reduce. In particular, cutoff effects may be larger in the 
presence of sea-quarks [8]. With Wilson fermions, even in their 0(a) improved version 
that we employ, it is well known that the chiral limit does not commute with the con- 
tinuum limit, implying that at fixed lattice spacing a cutoff effects become large below 
some quark mass. It is therefore important to control cutoff effects as one proceeds to 
simulate deeper in the chiral regime. 

In the quenched work [9j, rather accurate results were obtained in the pseudoscalar 
and vector channels using the Schrodinger functional. In this paper we carry over this 
computational setup to the Nf = 2 theory. The accuracy achieved [9] on masses was 
comparable to the calculations performed with periodic boundary conditions, and for 
decay constants the Schrodinger functional even proved to be the superior method. This 
is different when dynamical fermions are present. As shown in [TO] multi-pion excited 
states contribute significantly. For a computation of ground state masses and matrix 
elements they have to be supressed by a rather large time extent of the Schrodinger 
functional - in particular when the quark mass is low. In this situation it is more 



practical to employ (anti)periodic boundary conditions with the associated translation 
invariance in time. We can nonetheless use our simulation results to perform a first 
scaling test of the N{ = 2 0(a)-improved theory at low energies. Note that at high 
energies and correspondingly small lattice spacings excellent scaling has been seen [IH 
I12j . Besides the scaling test we give some details of our simulations including the 
algorithmic performance (section 2). 

2 Lattice simulations 

Our discretization consists of the Wilson gauge action and the non-perturbatively 0(a) 
improved Wilson quark action, with Cgw given in [13j. The algorithm and solver used in 
the present simulations have been described in some detail in [l3l|l5]. Using the notation 
of [16] for the hopping terms of the Dirac operatoi|l|, we recall the Schur complements 
of the hermitian Dirac operator with respect to asymmetric and symmetric even-odd 
preconditioning Qa, Q 

Q^ = c75(Moo-MoeM-iMeo), Q = M^-^Qa , c=(1 + 64k2)-i. (1) 

The action then reads 

S = Sg + Spf + 5det, (2) 



with 



'S'pf = (Ao 



QQt + pIm-,^] <Po + 4 [po^ + Qf] <t>i (3) 

5det = (-2)logdetMee + (-2) log dct Moo, (4) 



and Sg is the plaquette action. The determinants appearing in 5det are taken into 
account exactly. 

In Tab. [Hand Tab. [2] we list the simulations discussed in this paper. The reference 
length scale L* is defined through (f'{L*) = 5.5, where g is the Schrodinger functional 
coupling, and the values it assumes at the relevant bare couplings were presented in [lOj . 
For an estimate of L* in fermis, one may use the result a = 0.0784(10)fm at /5 = 5.3 [1], 
yielding L* ~ 0.6fm. 

Renormalization is carried out non-perturbatively in the SF at the scale p^en = 
l/Li-en, where 5^(Li.en) = 4.61. The values of the renormalization factor Zp of the pseu- 
doscalar density are taken from [T7| , while the values of the renormalization factor Za of 
the axial current differ from [T7| . They are presently re-evaluated using a Ward identity 
in a 1.8 fm x (1.2 fm)^ Schrodinger functional where the O(a^) effects are significantly 
smaller than before. In the table we list our preliminary numbers [18], which are not 
expected to change by more than the quoted errors. 

2.1 Stability and the spectral gap 

The spectral gap p of the Hermitian Dirac operator was used in [TDj as a tool to diagnose 
the stability of the HMC algorithm. We define 

p = - — -Tuvn \ VA A is an eigenvalue of QQ r > (5) 



^Moo, Mcc correspond to 1 + Too and 1 + Too respectively in [161115) . 



P {L/af X T/a 



L*/a 



Zk 



Zp 



Ax 


5.5 


323 X 42 


0.13630 


10.68(15) 


0.805(5) 


0.5008(70) 


Bi,B[ 






0.13550 








B2 


5.3 


243 X 32 


0.13590 


7.82(6) 


0.781(8) 


0.4939(34) 


B3 






0.13605 








B4 






0.13625 









C2 



5.2 



16^ X 32 
24^ X 32 



0.13568 
0.13568 



6.51(12) 0.769(12) 0.4788(5) 



Table 1: Simulation parameters. We use L*, defined by g'^{L*) = 5.5, as a reference scale. 
The renormalization factor of the axial current [171 I18j . Za, and of the pseudoscalar 
density [TT] at scale fi^en are listed. 



mol. dyn. JVrep ■ not 



Po 



(iv(°) 



CG/ 



(iV,« 



CG/ 



-Pa, 



Ai 


[LF; 2; 5; 50] 


1 


4340 


0.019803 


170 


824 


88% 


Bi 


[SW; 2; 1; 64] 


2 


2400 


0.0300 


100 


482 


91% 


B\ 


[SW; i; 1; 16] 
[SW; |; 1; 16] 


2 


1750 


0.0300 


100 


485 


90% 


B2 


2 


1900 


0.0300 


102 


729 


90% 


B3 


[LF; 2; 5; 50] 


2 


2600 


0.019803 


143 


905 


91% 


Bi 


[LF; 2; 5; 50] 


2 


1448 


0.0180 


155 


1195 


87% 


Ci 


[LF; 2; 5; 64] 


1 


6500 


0.0198 


179 


791 


96% 


C2 


[LF; 2; 5; 80] 


2 


2080 


0.0198 


184 


1086 


94% 



Table 2: Algorithmic parameters of the simulations. The molecular dynamics is char- 
acterized by [Integrator; r; (5Ti/(5To;T/(5ri], where the integrator can be 'leap-frog' or 
'Sexton- Weingarten' and subscripts refer to the two pseudofermions in use. For the 
gauge force, the SW integrator with (5ro/(5rg = 4 is used in all cases, and {Nqq) is 
the number of conjugate-gradient iterations used to solve the symmetrically even-odd 
preconditioned Dirac equation during the trajectory. 



03 




0.03 



a|x 



Figure 1: Histogram of ft for two different spatial volumes, simulations Ci and C2. The 
median is indicated in each case by the vertical dashed line. 



normalized such that it is given by the quark mass in the free theory with periodic 
boundary conditions. Since the only term that can potentially lead to unbounded fluc- 
tuations of the molecular dynamics forces is associated with Q, a sufficient condition 
for the stability of the algorithm is for the distribution of fi to be well separated from 
the origin. We remark that fi and /i (which was considered in [19]) cannot be directly 
compared on a quantitative level as they differ by the boundary conditions in the time 
direction and due to our (symmetric) even-odd preconditioning. We obtained fi by com- 
puting the lowest eigenvalue of QQ' using the algorithm of [20]. Figure (D displays the 
histogram of fi for simulations Ci^2- There is a clear separation of the median of the dis- 
tribution from the origin, but in a few cases in the course of the simulations eigenvalues 
as small as a third of this value were seen. 

We consider now the variance a^ of /i. In |19j . a measure a of the width of the fi 
distribution was found to approximately satisfy aay/l^Tja^ « constant. In the subset 
of our simulations where we computed fi, we find 



a 



VZ3T/C 



1.437(64) Ai 
1.268(23) Ci 
1.477(33) C2 



(6) 



varying only by about 15%. 



2.2 Autocorrelation times 

We compile observed integrated autocorrelation times Tint [22] in Tab.Ofor five quantities 
discussed and defined in detail in the next section. The dependence of the autocorre- 
lation times on the trajectory length was discussed previously [33]. Here we note that 
while there is a tendency for the autocorrelation time of the plaquette to decrease when 
the lattice spacing is increased, for other observables the opposite trend appears to be 
present underlining that autocorrelations have to be monitored for each observable sep- 
arately. The most important information in Tab. [3] is that all autocorrelations are small 
compared to the length of the runs (cf. Tab. [2]). Error estimates are hence trustworthy. 



Tint[0] 


P 


m{T/2) 


<siT/2) 


<AT/2) 


Feff(T/2) 


rnYAT/2) 


Geff(T/2) 


Al 


5.0(9) 


4.9(9) 


11(3) 


21(6) 


10(2) 


40(10) 


23(7) 


Bi 


13(3) 


5.5(9) 


7(1) 


16(4) 


4.2(7) 


23(7) 


11(3) 


B\ 


6(1) 


6(1) 


10(2) 


22(7) 


14(4) 


24(8) 


12(3) 


B2 


4.1(7) 


4.1(7) 


10(3) 


14(4) 


8(2) 


23(7) 


24(8) 


Bs 


9(2) 


3.9(6) 


4.7(7) 


11(2) 


6(1) 


11(3) 


11(2) 


Ba 


8(2) 


5(1) 


6(1) 


7(2) 


4.6(9) 


15(5) 


8(2) 


Cx 


9(2) 


5.3(8) 


5.2(8) 


5.1(8) 


4.7(7) 


4.9(7) 


5.6(9) 


C2 


11(3) 


6(1) 


6(1) 


7(1) 


3.9(6) 


6(1) 


6(1) 



Table 3: The integrated autocorrelation times for the plaquette, the current quark mass, 
the effective pseudoscalar mass and decay constant, and the effective vector mass. The 
unit is molecular dynamics time, i. e. trajectories times the length of the trajectory. For 
a precise definition of the observables see the following section. 



3 Scaling test 

In this section, which represents the central part of this paper, we investigate the cutoff 
effects on a number of non-perturbatively renormalized quantities. In order to keep 
systematic effects due to a varying volume negligible, we compare series of simulations 
in a fixed (but quite large) volume on a physical scale. More precisely we determine 
L/L* = 3.00(4), 3.07(3) and T/L* = 3.93(4), 4.09(3) on the A and B lattices. At 
P = 5.2, the volumes came out less uniform, L{Ci)/L* = 2.46(5), L{C2)/L* = 3.69(6) 
and T{Ci)/L* = 4.92(10). We shall discuss how to correct for these small mismatches 
after introducing the finite volume observables of this study. 

They are extracted from the zero spatial momentum boundary-to-bulk correlation 
functions, fA{xo), fp{xo) in the pseudoscalar channel, ky{xo) in the vector channel 
and the boundary-to-boundary pseudoscalar correlator /i [9j. We include the 0(a) 
improvement term proportional to ca 121] in /a,i = /a+oca^o/p. Effective masses 
and decay constants 

(7) 
(8) 
(9) 



m^fj(xo) = 


: -i(5o*+5o)log(/A,i(xo)) 


mfg(xo) = 


. -l{d*+do)log{fp{xo)) 


"irfr(^o) = 


. -i(a* + ao)iog(fcv(xo)) 





am 


am^ 


aml^ 


am\ 


aF^S 


a^G,s 




Za (l+bf^ama,) 


Zp (l+6pamq) 


Ai 


0.015519(37) 


0.1800(20) 


0.1793(15) 


0.2821(50) 


0.05999(42) 


0.0629(10) 


Bi 


0.03388(12) 


0.3272(18) 


0.3236(16) 


0.4520(35) 


0.09451(41) 


0.1507(14) 


B2 


0.019599(95) 


0.2391(35) 


0.2406(19) 


0.3953(51) 


0.08442(68) 


0.1267(22) 


B3 


0.01460(11) 


0.2118(24) 


0.2066(17) 


0.3647(35) 


0.07714(60) 


0.1170(13) 


Bi 


0.00727(14) 


0.1423(55) 


0.1528(20) 


0.3058(69) 


0.0698(11) 


0.0985(15) 


Ci 


0.01401(21) 


0.2173(55) 


0.2338(24) 


0.4354(60) 


0.0877(13) 


0.1637(25) 


C2 


0.01442(14) 


0.2328(39) 


0.2261(15) 


0.4152(42) 


0.08773(67) 


0.1614(15) 


Ci 


0.01431(19) 


0.2286(97) 


0.2282(63) 


0.410(14) 


0.08772(61) 


0.1620(17) 



Table 4: Simulation results for the effective quantities evaluated at xq = T/2. The bare 
current quark mass has been averaged over T/3 < xq < 2T/3. The last line gives the 
interpolation of Ci, C2, including the corrections described in the text. 



Feff(xo) = -IZk 



fAixo) (1 + 6Aamq)exp(m^g(a;o)(xo - T/2)) 



Geflf(xo) 



-2Za (1 + bAamq 
2Zp (1 + bpaniq) 



{hm^^{xo)L^f' 

{-^'^^^/'^ ,,, at xo = T/2 

(/imA^(T/2)L3)^/^ 

/p(xo) exp(m^jj(xo)(xo - T/2)) m^g{xo)^^^ 



(10) 



(/i L^) 



1/2 



, , ,/p(T/2) mP (T/2)i/2 

2Zp(l+fepamq) ^^^ \[ _f;,/ ' at xq 



T/2 



(11) 



(/i T3)V2 

are related to (T-dependent) masses and matrix elements, 

"T-rffl^^o) ~ Mps ^ rn^ff{xo) , "^^^(xo) ~ My , Teff(xo) ~ Tps , Ges{xo) ~ Gps -(12) 

These relations hold in the limit of large xq and T up to correction terms [9] 

Oeff (xo) = + 7J0 exp(-(^i - Mps) xo) + fio exp(-S2 (T - xq)) + . . . , (13) 

where the coefficients r]o and fio sue ratios of matrix elements, Ei is the energy of the 
first excitation in the zero momentum pion channel and E2 in the vacuum channel. For 
not too small L and not too large Mps we expect Ei ^ 3 Mps and E2 ^ 2 Mps . Our 
results for the effective observables at xq = T/2 are listed in Tab. S] together with the 
bare current quark mass m stabilized by averaging over T/3 < xq < 2T/3, 



1 



n2 



m 



n2- ni + 



- ^ wi(xo) , ni > T/3a , ^2 < 2T/3a 



m(xo) 



xo/a=ni 

i(9f5 + (9o)/a(xo) + CA a9^9o/p(xo) 



2/p(xo) 



(14) 



(15) 



The results at /3 = 5.3 can be compared directly to those of [3], shown in Tab. [5l 
for which the correction terms in Eq. (J13p can safely be neglected. In other words they 



6 





am 


a Mps 


aM\- 


aFps 


a2 Gps 




ZA(l+f>Aamq) 


Zp (l+6pamq) 


Di 


0.03386(11) 
0.01957(07) 
0.00761(07) 


0.3286(10) 
0.2461(09) 
0.1499(15) 


0.464(3) 
0.401(3) 
0.344(9) 


0.0949(13) 
0.0815(10) 
0.0689(13) 


0.1512(20) 
0.1260(16) 
0.1017(24) 



Table 5: Observables from fits of [3] i.e. xq^T — > cxd. Input parameters /?, k and L/a 
match those of lattices Bi,B2,B4^; note that D4 has been renamed here compared to [3]. 



correspond to xq, T — > cxd. This allows us to estimate the effects due to T{C) > T{A) ps 
T{B) in addition to those coming from the mismatch in L. 

1. For the matrix elements Fefr , Ggg no systematic differences between B and D 
lattices are visible. No correction due to T is necessary. We just interpolate the 
Ci and C2 results in L to L/L* = 3 using the Ansatz ai + 02 L~^/2e~*^Psi^ 
with Mps the pion mass on the larger volume. A small systematic error is added 
linearly to the statistical one. It is estimated by comparing with the result from 
an alternative interpolation with a'l + 02 L^^. 

2. We observe | m^g{B)/ m^g{D) — 1| < 0.03 without a systematic trend as a function 
of the quark mass. We take this into account as a systematic error of 2% on 'm^g{C) 
anqj subsequently we interpolate in L as in 1. The numbers for m^g are not used 
further. 

3. Finite T effects are not negligible in the vector mass {m^ff{B)/ m^g{D) — 1 ~ 
—0.10 ... — 0.03). We thus first perform a correction for the finite T effects using 
fits to Eq. (dSD with Ei = 2(Mps^ +(27r/L)2)V2^ E2 = 2Mps. A systematic error 
of 50% of this correction is included for the result. Next the finite L correction is 
performed as above. 

The interpolated values are included in Tab. H] as "simulation" Cj. After these small 
corrections we are ready to look at the lattice spacing dependence of our observables. To 
this end the necessary renormalization factors are attached (with perturbative values 
for 6a)&p [24J) and we form dimensionless combinations by multiplying with L* . At 
lowest order in the quark mass expansion (in large volume), one has Mpg oc m. It 
is thus natural to consider [rn^g L*]"^ /[rh{ij,ren)L*] instead of the quark mass itself. We 
choose fh renormalized non-perturbatively in the SF scheme at scale /Uron = l/-^rcn where 
5^(Lren) = 4.61 [TT]. The quantities considered are shown in Fig. [5] as a function of the 
dimensionless [m^gL*]^. At /? = 5.3 we have a few quark-mass points. As a reference, 
these are locally interpolated in [m^^L*]^ with a second order polynomial. For masses 
lighter than in simulation B2, the interpolation involves the lightest three masses and for 
heavier ones, it involves the heaviest three masses. The two-sigma bands (±2(t) of these 
interpolations are depicted as dotted vertical lines. Our results at the other /3- values 
are seen to be in agreement with these error bands, which are generally around 5%, but 



^From Eq. (|13|l this finite T effect scales with exp(— Mps T), yielding a reduction of 3% by a factor 
[1 — exp(— Mps L*)] when one considers the difference between T ^ 5L* and the target T = 4L* . 
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Figure 2: Dimensionless renormalized finite volume observables as a function of 



[<sL*?. 



From top to bottom Gcs{L*) 



m 



cff 



L* 



KffL*]V[7fi(^rcn)i1/15 



4FeffL* 

are shown. Squares, circles and triangle are for j3 = 5.2 , 5.3 , 5.5 respectively. Effective 
quantities are at xq = T/2. The dotted band is an interpolation of the 13 = 5.3 data as 
described in the text. 
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Figure 3: The effective pseudoscalar masses m^^ and m^^ in simulations B2 and Ai. Tlie 
liorizontal error bars are sliown on some of the points only for clarity. The horizontal 
line is to guide the eye. The vertical line indicates the middle of the B2 lattice. 



10% for 



m 



eflf 



L*\ /[m{^T.c-a)L*\ after all errors are included. Even if the precision is not 



very impressive, large cutoff effects are clearly absent. 

So far we have discussed the scaling of the ground state properties for a given symme- 
try channel. We now turn to the size of cutoff effects affecting excited state contributions 
to the correlators. Figure O compares the effective pseudoscalar masses m^j and rn^^ 



in simulation Ai and B2. The large size of the excited state contributions [10], while a 
drawback in extracting ground state properties, means that these functions are rather 
sensitive to the aforementioned cutoff effects. Because the Ai time extent is shorter by 
4(1)%, on this figure we have separately aligned the two boundaries of lattice Ai and 
i?2- We observe that the two data sets are consistent within uncertainties well before 
the function flattens off. With the exception of m^g for xq < T/2, the agreement sets 
in at a distance to the closest boundary of about L*, where it is easily seen that several 
excited states contribute significantly to the correlation functions. Altogether this figure 
is evidence that the masses and matrix elements of the first excited state in both the 
pion and vacuum channel have scaling violations not exceeding the few percent level. 
But higher states can have rather significant discretization errors. 

4 Conclusion and an updated value of A^ 

We carried out a finite size scaling test of the standard non-perturbatively 0(a)-improved 
Wilson theory with two flavors of dynamical fermions. In contrast to previous 



indications [8], cutoff effects are ratlier small in the present situation where the linear 
extent of the volume is around 1.6 fm. In fact within our precision of about 5% (collecting 
all errors) for effective masses and matrix elements, no a^ effects are visible. Continuum 
extrapolations of data from (say) 0.08 fm < a < 0.04 fm lattices which can nowadays be 
simulated [27\ I28j . seem very promising. Such a programme has been initiated [29j. A 
complementary effort |30) uses the twisted mass regularization of QCD [31j . Also in this 
case linear a-effects are absent [32] and the O(a^) effects appear to be moderate [55] . 

Finally we exploit the increased confidence in the scaling behavior of the simulated 
lattice theory to slightly refine our earlier estimate of the A-parameter. In |12] the 
product L* Ay^ = 0.801(56) was computed non-perturbatively in the two-flavor theory. 

Setting the scale through tq = 0.5 fm the value A^ = 245(16) (16) MeV was obtained 
emphasizing that more physical observables should be used in the future to set the scale. 
Given the quality of scaling observed in the previous section, it seems safe to assume 
that L*mK in the continuum limit differs by no more than 5% from its value at /3 = 5.3 
where rriY^a = 0.197(10) from [L,^ and L* /a = 7.82(6) [lOj are knowno We then obtain 

(2') f2) 

K^-L/rriK = 0.52(6) or A^ = 257(26) MeV, where a 5% uncertainty for a possible 

IVi o iVi o 

scaling violation has been added to the error (in quadrature). The new estimate is a bit 
higher than the previous one [12j . 
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